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Line spectra of 68 Taurus T Tauri stars were obtained with the Herschel-PACS (Photodetector Array Camera &■ Spectrometer) instrument as part of 
the GASPS (Gas Evolution in Protoplanetary Systems) survey of protoplanetary discs. A careful examination of the linescans centred on the [OI] 
63. 18 /jm fine-structure line unveiled a line at 63.32 /im in some of these spectra. We identify this line with the 8 is — * 7ov transition of ortho-water. 
It is detected confidently (i.e., > 3<x) in eight sources, i.e., ~ 24% of the sub-sample with gas-rich discs. Several statistical tests were used to search 
for correlations with other disc and stellar parameters such as line fluxes of [O I] 6300 A and 63. 18 yum; X-ray luminosity and continuum levels at 
63 /im and 850 yum. Correlations are found between the water line fluxes and the [O I] 63. 1 8 yum line luminosity, the dust continuum, and possibly 
with the stellar X-ray luminosity. This is the first time that this line of warm water vapour has been detected in protoplanetary discs. We discuss 
its origins, in particular whether it comes from the inner disc and/or disc surface or from shocks in outflows and jets. Our analysis favours a disc 
origin, with the observed water vapour line produced within 2-3AU from the central stars, where the gas temperature is of the order of 500-600K. 

Key words. Stars: formation, astrobiology, astrochemistry, line: identification, molecular data, protoplanetary discs. 
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1. Introduction 

Discs are natural by-products of star formation and the birth- 
places of planets. One of the key questions intimately linked 
with planet formation and the concept of planet habitability is 
how much vapour and icy water is present in discs and how it 
is radially distributed. However, it is only recently that observa- 
tions of water vapour in discs have become possible. 

ICarr & Naiital (I2008I) . using the Spitzer InfraRed 
Spectrograph (IRS), reported a rich molecular emission- 
line spectrum dominated by rotational transitions of hot water 
from the disc of AA Tau. They concluded that the molecular 
emission seen in the mid-IR has a most likely origin w ithin the 
2-3 AU inner regions of the disc. ISalvk et al.l (|2008|) detected 
water emission in the 10-20 yum region with Spitzer-IRS, as well 
as water and hydroxyl emission around 3 pm with NIRSPEC on 
Keck II, for both DR Tau and AS 205 A. The emission comes 
from the disc atmospheres of the objects and the excitation 
temperatures were found to be (~ 1 000 K), which is typical of 
terrestrial planet formation regions. IPontoppidan et al.l (12.0 1 Obh 
performed a survey for more protoplanetary discs in Ophiuchus, 
Lupus, and Chamaeleon, again with Spitzer-IRS, and concluded 
that the presence of mid-IR molecular emission lines, including 
those of water, is a common phenomenon in discs around 
Sun-like stars. Also, IPontoppidan et al.l (1201 Oal) presented a 
sample of ground-based observations of pure rotational lines 
of water vapour in the protoplanetary discs of AS 205 A and 
RNO 90 that was analysed to measure line widths of 30-60 km 
s , which is consistent with an origin in a disc in Keplerian 
rotation at a radius of ~ 1 AU, and gas temperatures in the range 
500-600 K. 

The Herschel Space Observatory dPilbratt et al.ll2010T) has 
opened the far-IR window with unprecedented sensitivity, al- 
lowing astronomers to survey the atomic and molecular gas 
content of disc regions that are cooler than those probed by 
Spitzer and ground-based instruments. Predictions of the water 



lines d etec table by Herschel ca n be found in Cernicharo et al.l 
(120091) and lWoitke etail (l2009bh . Irfogerheiide et all (1201 ll) pre^ 
sented the first detection of cold water in the disc of the ~10 
Myr old TW Hya, using Herschel-HIFI. The Herschel Open 
Time Ke y Program GASPS (G as Evolution in Protoplanetary 
Systems) dMathews et al.ll2010T) is conducting a survey to mea- 
sure gas lines and continuum in ~250 discs around low- 
and intermediate-mass stars with ages in the range l-30Myr 
with PACS, the Phot odetector Array Camera & Spectrometer 
dPoglitsch et al.ll2010l) . 

In this letter, we report the first detection of the 0-H2O line at 
63.32 /im in a subsample of protoplanetary discs around T Tauri 
stars in the 1-3 Myr old Taurus star forming region. 

2. Observations and data reduction 

This study is based on a sample of 68 classical and weak-line 
T Tauri stars from the Taurus star forming region with spectral 
measurements from Herschel-PACS centred at the wavelength 
of the [OI] 3 Pi — > 3 P2 63.184 pm line. The Taurus star-forming 
region is one of the main targets for the study of protoplanetary 
systems, because it is among the nearest star-forming regions 
(d = 140 pc) with a well-known populati on of more than 30 
young stars and brow n dw arfs according to Kenvon et al.l d2008l) . 
iLuhman et aD (1201 Oh . and lRebull etaiT(l2010l) . 

The observations described in this letter are part of the 
Herschel Open Time Key Programme GASPS [P.I. W. Dent], 
dMathews et al.ll2010l see), a flux-limited survey devoted to the 
study of the gas and dust in circumstellar systems around young 
stars. The survey focuses on the detection of the [OI] emission 
at the 63.18 pm feature. For this study, we analyzed 68 stars 
with spectral types ranging from late F-early G to mid M. The 
PACS spectral observations were made in chop/nod pointed line 
mode. The observing times ranged from 1215 to 6628 seconds, 
depending on the number of nod cycles. The data were reduced 
using HIPE 7.0.1751. A modified version of the PACS pipeline 
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Fig. 1. Spectra for the objects with a 63.32jum feature detection (> 3<x). The red lines indicate the rest wavelength of the [OI] and 
0-H2O emission. 



was used, which included: saturated and bad pixel removal, chop 
subtraction, relative spectral response-function correction, and 
flat fielding. Many observations suffer from systematic pointing 
errors, in some cases as large as 8", and are always shifted to the 
East. This is due to a plate scale error in the star tracker, which 
is normally negligible except in areas where the tracked stars 
are asymmetrically distributed within the field, as in Taurus. The 
mis-pointing translates into systematic small shifts in the line 
centre position. When the star was well centred within a single 
spaxel, we extracted the flux from that spaxel and applied the 
proper aperture correction. When the flux was spread over more 
than one spaxel, we co-added the spaxels. 



Table 1. Line positions and fluxes from PACS spectra. 



Name 


Sp. type 


Ah 2 o - \ai] 


[OI] flux 


o - H 2 flux 






fim 


10-" W/m 2 


10-" W/m 2 


AATau 


K7 


0.140 


2.2 + 0.13 


0.80 ± 0.13 


DLTau 


K7 


0.141 


2.4 + 0.15 


0.65 ±0.14 


FSTau 


M0 


0.139 


37 ± 0.26 


2.00 ± 0.33 


RYTau 


Kl 


0.139 


10 ± 0.42 


1.95 ± 0.38 


TTau 


K0 


0.138 


830 ± 0.75 


27.8 ± 0.7 


XZTau 


M2 


0.139 


32.2 ± 0.48 


2.11 ± 0.48 


HLTau 


K7 


0.142 


54.3 ±0.71 


8.14 ± 0.80 


UY Aur 


M0 


0.139 


33.6 ± 0.37 


1.90 ± 0.32 



Notes. All spectral types from the compilation ofLu hman et alj j201Cf) . 



3. Results and discussion. 

Among the sample of 68 Taurus targets studied in this letter, 
33 have discs that are rich in gas. These 33 all show the [OI] 
3 Pi — > 3 P2 line in emission at 63.18 fim (signal-to-noise ratio 
> 3, with values ranging from 3 to 375). In 8 of these 33 tar- 
gets (~ 24%), an additional fainter emission-line at 63.32 yum is 
detected (Fig.l). We computed 63.32 //m line fluxes by fitting 
a gaussian plus continuum curve to the spectrum using DIPSO 
Q. To improve the line fitting, the noisier edges of the spectral 
range were removed (i.e., A < 63.0 and A > 63.4). The re- 
sults are listed in Table [TJ where we report the peak position 
of the feature with respect to the observed wavelength of the 
[OI] 63.18 /jm line. According to these fits, the peak of the 
feature is at A = 63.32 /im. The FWHM is 0.020 jum, i. e., 
the instrumental FWHM for an unresolved line. We identify the 
feature as the ortho-FLO 8 ig — > 7 07 transition at 63.324 //m 
(Eupper Levei= 1070.7 K, Einstein A=1.751 s~ l ) since no other 
abundant species emit strongly at or close to the observed 
wavelength of the fe ature. This water feature was observed by 
iHerczee et all (1201 ll) in the outflow of NGC 1333 IRAS 4B. 
The o-FLO emission is only present in spectra with [O I] detec- 
tions. The targets FS Tau, HL Tau, and T Tau display extended 
emission in the [OI] 63.18 jt/m line, but only T Tau show hints 
of extended emission in the 63.32yum o-FLO line. T Tau is an 
exceptional object. It is a triple star system that drives at least 
two jets. It was t he most line rich PMS star observed by ISO 
dLorenzettil 120051) and resembles more a hot core than a PMS 
star, since the continuum emission is quite extended (~ 3.5 "at 

1 http://star-www.rl.ac.uk/docs/sun50.htx/sun50.html 
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Fig. 2. Plot of the 63.32 yum 0-H2O line luminosity versus [OI] 
63.18 /urn. Filled dots are detections, arrows are upper limits. 
Solid, dark grey arrows represent objects with non-detections 
spanning the same spectral range as the objects with detections. 
Light grey, empty arrows represent non-detections with other 
spectral types. 
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Fig. 3. Plot of the 63.32 //m ortho-H^O line flux versus 63yum 
continuum flux. Symbols as in Fig. [2] 



70 yt/m). It is impossible to tell how much of the line emission 
comes from the discs, from the outflows, or even from the sur- 
rounding envelope. 

To help us understand the origin of the 0-H2O emission, we 
compared the line intensity with several star and disc (jet) pa- 
rameters. We com puted survival analys is ranked statistics using 
the ASURV code ffeigelson & Nelsortfl98Msobe et al.lfl986h . 
The result of this analysis is summarised in Table [2] We also 
created random populations to test the validity of the results. 

The survival analysis shows a correlation between the 0-H2O 
line fluxes and the [O I] line fluxes at a significance level of 0.99 
(see Fig. 0. This relationship suggests that both lines have a 
similar origin. The H2O emission is correlated with the contin- 
uum emission at 63 //m, but at a significance level of 0.95 in the 
Spearman statistics (Fig. [3]). The 850 fim continuum flux can be 
used as a proxy for the amount of dust present in the disc. A 
survival analysis shows a possible correlation with the 850 yimi 
continuum flux, at a significance level of 0.95 in both Spearman 
and Kendall statistics, although a large scatter is present. A cor- 
relation with neither the stellar luminosity nor the spectral type 
is found. There seems to be a weak correlation with the slope 




10" 
Lx (erg/s) 

Fig. 4. H2O luminosity versus X-ray luminosity. Symbols as in 
Fig. [2] 

Table 2. Probabilities for correlations between 0-H2O line inten- 
sity and stellar/disc parameters. 



Observable 


Points n 


Spearman's prob. 


Kendall's prob. 


L[orj 


68 


0.0009 (0.4090 ) 


0.0000 (0.196) 


63 yum flux 


64 


0.0147 (0.2635 ) 


0.0002 (0.4567) 


850 fim flux (2) 


57 


0.0145 (0.2596) 


0.0131 (0.8496) 


T <3 » 


65 


0.1789 (0.1130) 


0.1980 (0.5535) 


Lx <3) 


65 


0.0225 (0.9912) 


0.0087 (0.6012) 


a(2yum - 8yum) <4) 


62 


0.026 (0.1547) 


0.0023 (0.3850) 


L[OI]6300A 


27 


0.1291 (0.3255) 


0.1008 (0.3450) 



Notes. The values obtained for random populations are shown in brack- 
ets. Accurate only if N > 30. (2): 850 /jm continuum fluxes from 
lAndrews & Williams! d2005h. (3): Val ues from lGudel etafl d2007h . (4): 
SEP slope from lLuhman et alJd20iol) . 



of the SED measured between 2 /mi and 8 yum, used as a proxy 
for the presence of hot dust. The significance of this correla- 
tion is dominated by T Tau, which has the highest 0-H2O flux 
in the sample. There is likely no link between mass loss rate 
and Lh,o- The JO I] luminosity at 6 300 A is proportional to the 
mass loss rate (lHartigan et al.lll995l) . Although the sample is too 
small to test this relationship conclusively, we note that while 
L[oi]630oA s P ans f° ur orders of magnitude, the H2O luminosity 
spans only one order of magnitude. Finally, the survival anal- 
ysis statistics points to a possible relationship with the X-ray 
luminosity L x (Fig. 3). While the Spearman probability for the 
real sample is only one order of magnitude smaller than for the 
random sample, the Kendall probability is two orders of magni- 
tude smaller. Inspection of Fig.|4]also shows that the 0-H2O line 
flux is detected only for sources with X-ray luminosities higher 
than 10 30 erg s _1 , which is consistent with photochemical disc 
models that show that far-IR line fluxe s increase signific antly 
above this X-ray luminosity threshold (lAresu et al ] 12011b . We 
note that log L x = 10 30 erg s _1 is above the median (mean) X-ray 
luminosity in Taur us (logL t = 29.8 (29.75) erg s _1 respectively, 
iGiidel et aD (120071) ). Interestingly, more than half of the sources 
with L x > 10 30 erg s do not display the 0-H2O line. This be- 
haviour may stem from either (1) the different shape of the X- 
ray spectrum (hardness ratio), (2) the duty cycle of the flares 
responsible for the high levels of X-ray fluxes, and/or (3) that X- 
rays are not the only driver of H2O chemistry and excitation, let 
alone any inner disc geometry and radiative transfer considera- 
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tions. We caution that the correlation with X-rays is significantly 
weaker when TTau is removed from the analysis. 

All the stars detected in 0-H2O are outflow/jet sources, al- 
though three of them AA Tau, DL Tau, and RY Tau do not show 
any excess in [OI], i.e. all the emission is consistent with coming 
from the disc (Howard et al., in prep). Two of these (AA Tau and 
DL Tau) are classified as outflow sources based on blue-shifted 
forbidden optical emission lines, although the emission is only 
slightly blue-shifted for AA Tau. However, in the sample there 
are also some prominent jet sources (£[oi]6300A > ^ 1 2 L©) that 
show no hint of emission from this 0-H2O line, such as DG Tau, 
where other H2O lines have been detected. However, these H2O 
lines are believed to originate in the outflow. Furthermore, the 
[OI] at 63.18 yt/mline is sometimes extended (Podio et al, in prep, 
Herczeg et al. 201 1), while we find the 0-H2O to be unresolved, 
suggesting a more compact origin for the 0-H2O line. 

We therefore assume that the 0-H2O emission at 63.32 fim 
originates from the disc and ask whether it comes from the 
same gas reservoir a s the hot H2O lines observed by Spitzer. 
ICarr&Naiita1(l2011h detected six out of eleven stars: AA Tau, 
BP Tau, DK Tau, GI Tau, RW Aur, and UY Aur. Our sample 
contains all of their detected sources. We det ected the 63.32 ^m 
H?0 em ission only in AA Tau and UY Aur. IPontoppidan et alJ 
d2010bl) reported H2O detections toward three out of eight stars 
in their sample: P R Tau, AA Tau, an d IQ Tau. We did no t detec t 
o- H^O in IO Tau. ISalvk et alJ d2008h . IPontoppidan et alJ d2009h 
and lMeiierink et al.ld2009l) argued that the Spitzer hot H2O emis- 
sion comes from the 0.1 to 1.0 AU annular region of the disc. 

Assumin g the same temperatu re, column density, and emit- 
ting areas as ICarr & Naiital d201 lb . i.e., 1AU, we computed the 
63.32//m LTE line flux in AA Tau and UY Aur. The model line 
emission is too optically thick to derive the H2O mass reliably. 
The size of the emitting region is instead estimated. In both 
cases, the model flux is ten times lower than measured. To re- 
cover the measured fluxes, the radius of the emitting area has to 
be about three times larger than the value they quote, i. e., 3.0 AU 
for AA Tau and 3.5 AU for UY Aur, a result that is consistent 
with the lower excitation temperature of the line we observed. 

A radiation thermo-chemical model of a typical T Tauri 
disc obtained with the ProDiMo code dWoitke et alJ l2009at 
iKamp et all 120101: iThi et all 1201 Ot lAresu et alJ 1201 lh predicts 
that the emission region of the 63.32 fim 0-H2O line is of the 
order of 3 AU, about five times larger than the Spitzer emission- 
line region (see Fig. IB.U . One particular 0-H2O line at 15.738 
fim was selected as a representative Spitzer mid-IR H2O emis- 
sion line. This model does not intend to fit any particular ob- 
ject. A large non-LTE H2O ro-vibrational model calculation 
was included to consistently calculate the Spitzer as well as 
the Hersch el H2O emission l ines by applying escape prob abil- 
ity theory ( Fa ure et al 1 120041 120071 iFaure & Josseiir3l2008l) . At 
~ 3 AU, the gas densities are high enough to excite the 63.32 
fim 0-H2O line and the dust temperature close to the mid-plane 
is low enough for H2O to freeze onto grains. 

The 63.32 fim line could provide the missing link between 
the mid-IR Spitzer detections of hot H2O vapour in T Tauri discs 
and the cold far-IR H2O lines observed with Herschel. Searches 
for the lower excitation H2O lines have proven to be less success- 
ful. To date, the only clear detection of cold H 2O from a disc is 
TW Hya using HIFI dHogerheiide et al.ll201 lh , and Kamp et al. 
(in prep.), using PACS. Possible reasons for the lower detection 
rate of cold H2O might be the much lower H2O abundances in 
the outer disc caused by the free ze-out of H7O and/o r significant 
(vertical) settling of icy grains dBergin et alJl2010T) . Observing 
the same molecule in transitions with very different excitation 



temperatures may trace it through a broader range of different 
radial zones in protoplanetary discs. 

4. Conclusions 

We have detected 0-H2O emission at 63.32 fim in 8 T Tauri stars 
in a sub-sample of 68 stars located in Taurus. The detection rate 
is ~24% in the sub-sample with gas-rich discs. The H2O emis- 
sion appears to be correlated with the continuum luminosities, 
the [OI] 63.18 fim line fluxes, and the X-ray luminosities. The 
gas temperature (500-600K) and density needed to excite the 
observed 0-H2O line suggest that the line is coming from the 
inner parts of the discs and from the upper layers of its atmo- 
sphere, where the disc is directly illuminated. The correlation 
with X-rays flux and the role of X-ray emission in heating the 
gas, in particular during flares, needs to be investigated further. 
The typical size of the emitting region is estimated to be r ~ 
3 AU, which is consistent with the typical location of the snow 
line in these objects. More effort is needed to detect several H2O 
lines simultaneously in more objects to understand the full radial 
distribution of H2O vapour in planet-forming discs. 
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Appendix A: Non detections 



GASPS is a flux-limited survey. All sources have detection lim- 
its for the 0-H2O line at 63.32 /urn around 10~ 17 W/m 2 , with an 
uncertainty of a factor of two both ways. The list of non detec- 
tions is given in table lA.ll 



Appendix B: Location of emission for H 2 lines 
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Fig. B.l. Location of the H2O emission for a typical disc model 
(M«=O.8M , L»=0.7 L Q , r eff =4400K, UV excess / uv =0.01, 
X-ray luminosity L x =10 30 erg s" 1 , fl in =0.1AU, fl O ut=300 AU, 
M disc =0.01 M©, dust/gas=0.01, e=-l, scale height H =l AU 
at ro=10 AU, flaring f3 = 1.1, astronomical silicate with uni- 
form dust size distribution a m ; n =0.05 /an, a max =l mm, and in- 
dex=3.5.). The 15.738 yum line is used as a typical HhO line de- 
tected by Spitzer. We have marked the radii where the cumula- 
tive line flux reaches 15% and 85% with vertical lines, for the 
four selected H2O lines. In addition, the horizontal lines mark 
the heights above the midplane where 15% and 85% of the line 
flux, from every vertical column, originate in. The encircled re- 
gions are hence responsible for 70% x 70% = 49% of the total 
line flux (for a pole-on disc). 
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Table A.l. Source list, spectral types, and 63.32 fim 0-H2O line 
fluxes. 



Name 


SP. Type 


- H 2 flux 


- 


- 


10"" W/m 2 


Anon 1 


MO 


<1.17 


BPTau 


K7 


<0.94 


CIDA2 


M5.5 


<0.89 


CITau 


K7 


<0.92 


CoKu Tau/4 


K3 


<1.37 


CWTau 


K3 


<1.26 


CXTau 


M2.5 


<0.94 


CYTau 


M1.5 


<1.41 


DETau 


Ml 


<1.23 


DFTau 


M2 


<1.28 


DGTau 


K6? 


<1.32 


DG Tau B 


<K6 


<1.25 


DHTau 


Ml 


<1.07 


DK Tau 


K6 


<0.53 


DL Tau 


K7 


<1.14 


DM Tau 


Ml 


<0.89 


DN Tau 


MO 


<0.77 


DO Tau 


MO 


<1.21 


DP Tau 


M0.5 


<0.92 


DQ Tau 


MO 


<0.98 


DS Tau 


K5 


<0.58 


FF Tau 


K7 


<1.11 


FM Tau 


MO 


<1.24 


FO Tau 


M3.5 


<0.96 


FQ Tau 


M3 


<0.91 


FT Tau 


- 


<1.25 


FW Tau 


M5.5 


<1.08 


FX Tau 


Ml 


<1.82 


GG Tau 


M5.5 


<1.24 


GH Tau 


M2 


<0.84 


GIK Tau 


K7 


<0.74 


GM Aur 


K7 


<1.21 


GO Tau 


MO 


<1.23 


Haro6-13 


MO 


<1.43 


Haro6-37 


Ml 


<1.84 


HBC 347 


Kl 1 


<1.69 


HBC 356 


K2 1 


<0.94 


HBC 358 


M3.5 


<1.29 


HK Tau 


M0.5 


<0.90 


HN Tau 


K5 


<0.64 


HO Tau 


M0.5 


<1.46 


HV Tau 


Ml 


<0.72 


IP Tau 


MO 


<0.83 


IQ Tau 


M0.5 


<1.38 


04158+2805 


M5.25 


<1.06 


LkCa 1 


M4 


<1.06 


LkCa3 


Ml 


<0.94 


LkCa4 


K7 


<0.98 


LkCa5 


M2 


<1.16 


LkCa7 


MO 


<0.88 


LkCa 15 


K5 


<0.52 


SAO 76428 


F8? 


<1.58 


SUAur 


G2 


<0.72 


UX Tau 


K5 


<0.99 


UZTau 


M2 


<0.60 


V710Tau 


M0.5 


<1.17 


V773 Tau 


K3 


<0.77 


V819 Tau 


K7 


<1.21 


V927 Tau 


M4.75 


<0.88 


VY Tau 


MO 


<1.30 



Notes. S pectral types from th e compilation by iLuhman et alj j201fj) . 
(1): from lHerbig & Belli d!988h 
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